Shock wave propagation in vibrofluidized granular materials

Shock wave formation and propagation in two-dimensional granular materials under vertical vibration are studied by digital high speed photography. The steepen density and temperature wave fronts form near the plate as granular layer collides with vibrating plate and propagate upward through the layer. The temperature front is always in the transition region between the upward and downward granular flows. The effects of driving parameters and particle number on the shock are also explored.Comment: 9 pages, 4 figures, submitted to PR

Mach cone in a shallow granular layer

We study the V-shaped wake (Mach cone) formed by a cylindrical rod moving through a thin, vertically vibrated granular layer. The wake, analogous to a shock (hydraulic jump) in shallow water, appears for rod velocities v_R greater than a critical velocity c. We measure the half-angle, theta, of the wake as a function of v_R and layer depth, h. The angle satisfies the Mach relation, sin(theta)=c/v_R, where c=sqrt(gh), even for h as small as one particle diameter.Comment: 4 pages 4 figure

Power-law velocity distributions in granular gases

We report a general class of steady and transient states of granular gases. We find that the kinetic theory of inelastic gases admits stationary solutions with a power-law velocity distribution, f(v) ~ v^(-sigma). The exponent sigma is found analytically and depends on the spatial dimension, the degree of inelasticity, and the homogeneity degree of the collision rate. Driven steady-states, with the same power-law tail and a cut-off can be maintained by injecting energy at a large velocity scale, which then cascades to smaller velocities where it is dissipated. Associated with these steady-states are freely cooling time-dependent states for which the cut-off decreases and the velocity distribution is self-similar.Comment: 11 pages, 9 figure

Persistent holes in a fluid

We observe stable holes in a vertically oscillated 0.5 cm deep aqueous suspension of cornstarch for accelerations a above 10g. Holes appear only if a finite perturbation is applied to the layer. Holes are circular and approximately 0.5 cm wide, and can persist for more than 10^5 cycles. Above a = 17g the rim of the hole becomes unstable producing finger-like protrusions or hole division. At higher acceleration, the hole delocalizes, growing to cover the entire surface with erratic undulations. We find similar behavior in an aqueous suspension of glass microspheres.Comment: 4 pages, 6 figure

Onset of Patterns in an Ocillated Granular Layer: Continuum and Molecular Dynamics Simulations

We study the onset of patterns in vertically oscillated layers of frictionless dissipative particles. Using both numerical solutions of continuum equations to Navier-Stokes order and molecular dynamics (MD) simulations, we find that standing waves form stripe patterns above a critical acceleration of the cell. Changing the frequency of oscillation of the cell changes the wavelength of the resulting pattern; MD and continuum simulations both yield wavelengths in accord with previous experimental results. The value of the critical acceleration for ordered standing waves is approximately 10% higher in molecular dynamics simulations than in the continuum simulations, and the amplitude of the waves differs significantly between the models. The delay in the onset of order in molecular dynamics simulations and the amplitude of noise below this onset are consistent with the presence of fluctuations which are absent in the continuum theory. The strength of the noise obtained by fit to Swift-Hohenberg theory is orders of magnitude larger than the thermal noise in fluid convection experiments, and is comparable to the noise found in experiments with oscillated granular layers and in recent fluid experiments on fluids near the critical point. Good agreement is found between the mean field value of onset from the Swift-Hohenberg fit and the onset in continuum simulations. Patterns are compared in cells oscillated at two different frequencies in MD; the layer with larger wavelength patterns has less noise than the layer with smaller wavelength patterns.Comment: Published in Physical Review

Enskog Theory for Polydisperse Granular Mixtures. I. Navier-Stokes order Transport

A hydrodynamic description for an $s$-component mixture of inelastic, smooth hard disks (two dimensions) or spheres (three dimensions) is derived based on the revised Enskog theory for the single-particle velocity distribution functions. In this first portion of the two-part series, the macroscopic balance equations for mass, momentum, and energy are derived. Constitutive equations are calculated from exact expressions for the fluxes by a Chapman-Enskog expansion carried out to first order in spatial gradients, thereby resulting in a Navier-Stokes order theory. Within this context of small gradients, the theory is applicable to a wide range of restitution coefficients and densities. The resulting integral-differential equations for the zeroth- and first-order approximations of the distribution functions are given in exact form. An approximate solution to these equations is required for practical purposes in order to cast the constitutive quantities as algebraic functions of the macroscopic variables; this task is described in the companion paper.Comment: 36 pages, to be published in Phys. Rev.

Scaling, Multiscaling, and Nontrivial Exponents in Inelastic Collision Processes

We investigate velocity statistics of homogeneous inelastic gases using the Boltzmann equation. Employing an approximate uniform collision rate, we obtain analytic results valid in arbitrary dimension. In the freely evolving case, the velocity distribution is characterized by an algebraic large velocity tail, P(v,t) ~ v^{-sigma}. The exponent sigma(d,epsilon), a nontrivial root of an integral equation, varies continuously with the spatial dimension, d, and the dissipation coefficient, epsilon. Although the velocity distribution follows a scaling form, its moments exhibit multiscaling asymptotic behavior. Furthermore, the velocity autocorrelation function decays algebraically with time, A(t)= ~ t^{-alpha}, with a non-universal dissipation-dependent exponent alpha=1/epsilon. In the forced case, the steady state Fourier transform is obtained via a cumulant expansion. Even in this case, velocity correlations develop and the velocity distribution is non-Maxwellian.Comment: 10 pages, 3 figure

Patterns and Collective Behavior in Granular Media: Theoretical Concepts

Granular materials are ubiquitous in our daily lives. While they have been a subject of intensive engineering research for centuries, in the last decade granular matter attracted significant attention of physicists. Yet despite a major efforts by many groups, the theoretical description of granular systems remains largely a plethora of different, often contradicting concepts and approaches. Authors give an overview of various theoretical models emerged in the physics of granular matter, with the focus on the onset of collective behavior and pattern formation. Their aim is two-fold: to identify general principles common for granular systems and other complex non-equilibrium systems, and to elucidate important distinctions between collective behavior in granular and continuum pattern-forming systems.Comment: Submitted to Reviews of Modern Physics. Full text with figures (2Mb pdf) avaliable at http://mti.msd.anl.gov/AransonTsimringReview/aranson_tsimring.pdf Community responce is appreciated. Comments/suggestions send to [email protected]

Diffusion in a Granular Fluid - Theory

Many important properties of granular fluids can be represented by a system of hard spheres with inelastic collisions. Traditional methods of nonequilibrium statistical mechanics are effective for analysis and description of the inelastic case as well. This is illustrated here for diffusion of an impurity particle in a fluid undergoing homogeneous cooling. An appropriate scaling of the Liouville equation is described such that the homogeneous cooling ensemble and associated time correlation functions map to those of a stationary state. In this form the familiar methods of linear response can be applied, leading to Green - Kubo and Einstein representations of diffusion in terms of the velocity and mean square displacement correlation functions. These correlation functions are evaluated approximately using a cumulant expansion and from kinetic theory, providing the diffusion coefficient as a function of the density and the restitution coefficients. Comparisons with results from molecular dynamics simulation are given in the following companion paper